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1 Objective 


The objective of this contract is to use the detailed particle energy spectra measured by the Los 
Alamos Magnetospheric Plasma Analyzer (MPA) on several satellites in geosynchronous orbit to 
model the potential of the chassis due to surface charging of those satellites, as derived also from 
MPA measurements. In particular, we seek to identify a relatively simple characterization of the 
full particle distributions that yields an accurate prediction of the observed charging under a wide 
variety of conditions. The three-dimensional computer codes that can compute spacecraft 
charging in a tenuous plasma (NASCAP/GEO [Katz, et al, 1979], SEE Spacecraft Charging 
Handbook [Katz, et al, 2000], and NASCAP-2K [ Davis et al, 2002]) assume that the 
environment can be adequately modeled by a Maxwellian or Double Maxwellian distribution 
function. We seek to verify this or determine an appropriate alternative formulation. 

2 Summary of Work Completed 

In preparation for this effort, we reviewed three previous closely related studies: [Garrett, et al., 
1980], [Lai, et al., 2001], and an internal study done by LANL. 

The measurements used in this study were taken by the LANL MPA, which is described in 
[Bame, et al., 1993], The determination of the spacecraft potential, the moments of the 
distribution, and the distribution function from the measurements are described in the following 
publications: [Lawrence, et al., 1999], [McComas, et al., 1993], [Thomsen, et al., 1996], 
[Thomsen, et al., 1999], Very simplified descriptions of some of this processing are included 
below in order to explain how the measurements are being used in this project. 

Seven datasets were constructed for this investigation. Each dataset record consists of an electron 
and an ion spin-averaged, differential flux spectrum, along with a time tag, the chassis potential 
derived from the ion line, and various moments of the distribution. The first two datasets consist 
of spectra from September 2001 for satellite 1994-084. These datasets were filtered for data 
quality and existence of an ion line. The first dataset was for spectra taken during eclipse while 
the second included times before and after eclipse. These datasets were used for construction and 
testing of the software needed to handle a large number of spectra at once. 

After the preliminary investigation was complete and the bulk of the software built, five 
additional datasets were constructed. All the results reported here use these five datasets. Four of 
the datasets are from the satellite 1994-084. This includes two months near the peak of the solar 
cycle, March 2001 and September 2001, and two months near the minimum, March 1996 and 
September 1996. The fifth dataset is from September 2001 for the LANL-97a spacecraft. This 
final dataset allows us to confirm that our results are valid for a different spacecraft of the same 
design. 

We developed software to examine a large number of spectra at once. The software, written in 
Java, does the following: 

• Read in a dataset, 

• Determine if the spacecraft is in eclipse or sunlight, 
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• Filter the dataset for sunlit/eclipse, presence of a potential computed from an ion line, 
and data quality, 

• Compute the moments, fluxes, and derived potential described in this report from 
measured fluxes and analytic fits to measured fluxes, 

• Fit a measured distribution function with an analytic form, and 

• Compute statistical measures. 

We created Excel workbooks to display the results as shown in this report and compute 
additional statistical measures. 

For the various datasets, we computed electron and ion density and temperature moments for the 
measured spectra. The results verify that the calculations done at SAIC using the Java software 
using the spin-averaged spectrum give the same results as those done at LANL using the full 
angular distribution, with the exception of the density of ions with measured energies over 
100 eV. The density and temperature moments of electrons with incident energies above 30 eV, 
the temperature moments of ions with incident energies above 100 eV, and the density moments 
of the ions with incident energies below 100 eV all agree well. The SAIC computed density 
moments of ions with incident energies above 100 eV are consistently lower than the LANL 
computed values. In addition to providing verification that the calculations are done correctly, 
doing these calculations allowed us to refine the calculational technique and become familiar 
with the data. 

For the various datasets, we computed incident electron and ion fluxes, secondary electron fluxes 
from incident electron and ion impact, backscattered electron fluxes, and the resulting net fluxes. 
We investigated how material properties, some choices in techniques to do the energy integral, 
and assumptions regarding the fraction of low energy particles escaping effect the calculated net 
current. 

We also investigated techniques to compute the photoemission fluxes. Geometric effects are 
important enough in the sunlit charging of the LANL spacecraft that current balance results in 
sunlight are more suggestive than quantitative. 

In order to gain a better measure of the quality of our flux computations, we computed the 
spacecraft chassis potential that we expect given incident electron and ion spectra for the various 
datasets. This approach provides a better prediction of chassis potentials than the measured 
temperature or the integrated measured flux between 9 and 47 keV. 

Using the measured spectra, we detemiined that computing the fluxes using a set of material 
properties for graphite carbon and a low-energy secondary electron escape fraction of 81%, gives 
computed potentials consistent with measured values. The estimated potential is within the a 
factor of 1 .5 of the measured potential 65% of the time and within a factor of 3 of the measured 
potential 87% of the time. 

Our goal is to provide guidance to those doing spacecraft surface charging calculations. We 
compared the results of calculations using various simple analytic fits with calculations using the 
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data. We developed an automated fitting procedure and computed moments, fluxes, and 
predicted potentials using Maxwellian, Double Maxwellian, and Kappa distribution functional 
forms. Christon, et al [1989] have shown that the Kappa distribution, which is a Maxwellian at 
low energies and a power law at high energies, often provides a good fit to the quiescent plasma 
sheet environment. (The spectrum is typically much more complex during disturbed conditions 
[Christon, et al, 1991].) 

The potentials computed using Maxwellian fits for both ions and electrons are within a factor of 
1 .5 of the measured potential 30 to 58% of the time and within a factor of 3 of the measured 
potential 61 to 71 % of the time, depending on the fitting procedure. Using a Kappa distribution 
for the electrons and a Maxwellian distribution for the ions gave results within a factor of 1.5 of 
the measured potential 65% of the time and within a factor of 3 of the measured potential 84% of 
the time. Potential predictions using the Kappa-Maxwellian fit combination gives results similar 
to those obtained from using the measured fluxes directly. Energy bins with spurious 
measurements must be discarded from the flux integral when using the measured fluxes directly. 
These energy bins and the energy range above the maximum measured can be included when 
using fit spectra. Using all energy bins from 1 eV to 100 keV does not improve the agreement. 

2. 1 Reports , Meetings, Presentation, and Publication 

This report includes the conclusions of the year-long study, as well as, all the techniques used 
and the results previously reported. Appendix A includes the quarterly reports for the LANL 
portion of this effort, including the final report. 

The kick-off presentation for this contract was held at SA1C in San Diego on January 29, 2002. 
The presentation summarized the intended work plan. 

Dr. Davis (Principal investigator) and Dr. Thomsen (Co-Investigator) were able to meet in 
person three times to discuss this contract. Dr. Thomsen traveled to SA1C in San Diego, CA, to 
attend the contract kick-off meeting and discuss the initial work on this contract. Dr. Davis 
traveled to LANL in Las Alamos, NM, where we had a productive full day meeting during which 
we reviewed how the data is handled, how the various integrals are calculated, and the meaning 
of the results. We also took advantage of an unrelated meeting that Dr. Thomsen had in Los 
Angeles, CA, to meet in San Diego, CA, for a half-day, where we discussed the appropriate high 
energy cut off for the integrals and improved the explanation of the justification for the technique 
used to compute the potential given the incident fluxes. 

A presentation at the Fall Meeting of the American Geophysical Union is included as Appendix 
B. 

A paper summarizing the results of this study will be submitted to the Journal of Geophysical 
Research. This final report and the paper will make the results of this study available to other 
researchers. 
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2.2 Conclusions 


The LANL dataset has proved to be a powerful tool for the investigation of spacecraft surface 
charging. The spectra provide the resolution and accuracy adequate for “post-diction” spacecraft 
charging calculations. 

It is necessary to include all the current components— incident electrons and ions, secondary 
electrons, backscattered electrons, and photoelectrons— to accurately “post-diet” chassis 
potentials from measured spectra. Eclipse calculations provide good accuracy. Calculations in 
sunlit conditions were significantly less successful. Estimates of photoemission and fraction of 
low energy electrons escaping are strongly dependent on geometric effects and sun angle. The 
quality of the data is high enough that with additional effort, we believe that the photoemission 
as a function of sun angle and chassis potential could be determined. This procedure would only 
work if the spacecraft surfaces are primarily conducting, as are most geosynchronous spacecraft. 
Potential “post-dictions” using a Kappa distribution to fit the incident electron spectrum and a 
Maxwellian distribution to fit the incident ion spectrum give results similar to “post-dictions” 
using the measured spectra directly. Better results would be obtained if additional intelligence 
were added to the ion flux fitting procedure. 

We recommend the addition of the Kappa distribution to the selection of options in the three- 
dimensional spacecraft charging codes. The addition of the capability of using measured fluxes 
would also be a useful addition to the NASCAP-2K software. 

3 Data 

The measurements used in this study were taken by Magneto spheric Plasma Analyzer (MPA) 
instruments built by the Los Alamos National Laboratory (LANL) and flown on a series of 
geosynchronous spacecraft. The MPA is a spherical-sector electrostatic analyzer with a bending 
angle of 60°. The field of view is divided into six separate but contiguous detectors covering the 
range of polar angle from about 25° to 155°. The satellite spin allows the instrument to view 360 
degrees in azimuth, which is divided into 24 sectors of 15 degrees. Thus in one spin the MPA 
views ~92% of the unit sphere, divided into six polar by 24 azimuthal view directions. The spin 
axis of the spacecraft points continuously at the center of the Earth, so the two polar angle 
detectors which view nearly perpendicular to the spin axis give very complete pitch angle 
coverage. While the spacecraft spins through a 15° azimuthal sector, the MPA plate voltage is 
swept through 40 logarithmically spaced energy channels ranging from ~40 keV/e down to ~1 
cV/e . A complete three-dimensional (40 energies x 24 azimuths x 6 polar angles) distribution is 
obtained in one 10-s spin. Since the same analyzer is used for ion and electron measurements (by 
changing the polarity of the plate voltage and channel electron multiplier bias [Bame et al., 
1993]), the ion and electron distributions are measured alternately. In 86 s, the instrument cycles 
through one three-dimensional electron distribution and two three-dimensional ion distributions, 
as described above, as well as three two-dimensional electron distributions and two high- 
angular-resolution modes [Bame et al., 1993]. 

The determination of the spacecraft potential, the moments of the distribution, and the 
distribution function from the measurements are described in the following publications: 
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[Lawrence, et al, 1999], McComas, et al., 1993], [Thomsen, et al., 1996], [Thomsen, et al., 
1999]. Simplified descriptions of some of this processing are included below in order to explain 
how the measurements are being used in this project. 

Seven datasets were constructed. Each dataset record consists of an electron and an ion spin- 
averaged, differential flux spectrum, along with a time tag, the chassis potential derived from the 
ion line, and various moments of the distribution. 

The first two datasets consist of spin-averaged, differential flux spectra from September 2001 for 
satellite 1994-084. The first dataset is confined to times when the spacecraft was between -0.6 
and +0.5 local time; for almost all the points the spacecraft was in eclipse. The second dataset 
consists of measurements taken when the spacecraft was between 22.4 local time to +2.0 local 
time; both sunlit and eclipse data are included. These initial datasets include only spectra for 
which a clear potential determination can be made from the ion line. These initial two datasets 
consist of 839 and 1 195 spectra respectively. These datasets were used for construction and 
testing of the software needed to compute moments and fluxes of a large number of spectra at 
once. 

After the preliminary investigation was complete and the bulk of the software built, five 
additional datasets were constructed. All the results reported here use these five datasets. Four of 
the datasets are from the satellite 1994-084. This includes two months near the peak of the solar 
cycle, March 2001 and September 2001, and two months near the minimum, March 1996 and 
September 1996. (Solar cycle 23 began in May 1996 and peaked in April 2000.) The fifth dataset 
is from September 2001 for the LANL-97a spacecraft. This final dataset allows us to confirm 
that our results are valid for a different spacecraft of the same design. 

The datasets include one record for each 86-second measurement period during the time period 
covered. The analysis software filters the data to keep only spectra for which the ion spectrum 
provides a distinct ion line to determine the potential along with sunlit/eclipse condition and 
acceptable data quality flags. Each record includes several derived quantities. The most 
important of these are the measured potential and its quality flag. Additionally, there are ion and 
electron densities and parallel and perpendicular temperatures, along with an estimated potential 
barrier height and its quality flag. The estimated potential barrier is the energy of the minimum 
electron count rate between 30 eV and 400 eV. Each record also includes a timestamp and the 
geocentric position from which the determination of sunlit or eclipse condition is made. Each of 
these quantities is further described in the section below on moments. 

3.1 Spectra Of A Charged Spacecraft 

Figure 1 and Figure 2 show the ion flux spectrum and the electron flux spectrum, respectively, 
for day 245.70276. (The fractional part of the day gives the time (GMT)). The points labeled 
“Flux at spacecraft” are the spin-averaged, measured fluxes. The points labeled “One-count” are 
the fluxes that would be measured if one particle hit one of the six detectors. This is an 
approximate value as the six detectors have different geometric factors and efficiencies. For this 
work, we treat flux values under twice the one-count rate as zero. 
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The energy channel in which the ion count rate increases dramatically determines the spacecraft 
potential. LANL computed values are used here. In Figure 1, the lowest energy channel with a 
significant count rate is 2450 to 3200 eV. The chassis potential is taken to be the geometric mean 
of the channel edges, -2800 V. An ion with nearly zero energy at infinity is accelerated to 2800 
eV by the time it reaches the detector, which is why no ions are seen below the 2450 to 3200 eV 
energy channel. The negatively charged spacecraft repels electrons. An electron with energy of 
2800 eV reaches the detector with zero energy, and lower energy electrons do not reach the 
detector at all. However, most of the electrons that reach the detector with energies of less than 
30 eV (and on highly charged spacecraft, sometimes even 200 eV) are generated on the 
spacecraft surface or in the spacecraft vicinity. The low energy peak seen in Figure 2 is electrons 
trapped by differential charging electric fields. 

Day 245.70276 with spacecraft at -2800 V 



Figure 1. Ion fluxes at day 245.70276, with spacecraft at -2800. 
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Day 245.70276 with spacecraft at -2800 V 



Electron energy (eV) 

Figure 2. Electron fluxes at day 245.70276, with spacecraft at-2800 V. The vertical line of triangles is locally 
created electrons. 


The accuracy of the energy of each measurement can be estimated as of the order of half the 
energy bin in which the measurement is made plus half the energy bin of the ion line. The 
accuracy of the flux is given by Poisson statistics to be the geometric mean of the measured spin- 
averaged flux and the one-count flux. 


3.2 Adjusting Spectra To Account For Chassis Potential 

The ambient fluxes can be derived from the fluxes at the spacecraft. The formula is derived from 
a consideration of the distribution function. The differential flux as a function of energy is related 

2Ee 2 

to the distribution function by the formula F(E) = — — f (E) , where F is the flux and f is the 

m 

distribution function. The measured energy is shifted from the energy at infinity by the potential; 
E m = E m + (f) , where the sign is for ions and the “+” sign is for electrons. 


Liouville’s theorem gives that “the density of system in the neighborhood of some given system 
in phase space remains constant in time.” [Goldstein, 1950a] Krall and Trivelpiece [1986] state 
Liouville’s theorem as “f(x, v, t) is constant along any particle trajectory.” For our purposes, we 
can state Liouville’s theorem as “if all possible particle trajectories that begin at the spacecraft 
end at infinity, the distribution function at the spacecraft is the same as the distribution function 
at infinity.” As we are treating the spacecraft as a uniform sphere, the potential varies only 
radially. All particles in a attractive radially symmetric potential field that varies more slowly 
than the inverse distance square connect to infinity. If the potential field varies faster than the 
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inverse distance squared, some particle trajectories both begin and end on the sphere due to 
angular momentum. [Goldstein, 1950b]. In a repulsive radially symmetric potential field, all 
particle trajectories connect to infinity. 


As the distribution function at the spacecraft is the same as the distribution function at infinity, 
we have that 


f-(E>) = f m (E m ) 
^(Ej = ^F m (E m ) 


F4E m ±9) = F m (E m ) 


f . \ 

l±± 

E J 


where the “+” sign is for ions and the sign is for electrons. 

Figure 1 and Figure 2 show the flux spectra shifted to infinity. The ion spectrum is fairly flat and 
the electron spectrum shows a spurious peak just above the spacecraft potential. The lower 
energy electrons were generated at the spacecraft surface and therefore do not reflect the 
properties of the spectrum at infinity. 

3.3 Spectra Of An Uncharged Spacecraft 

Figure 3 shows electron and ion flux spectra for an uncharged spacecraft (-5.4 V). The spectrum is 
definitely not Maxwellian. There is a low energy (about 7 eV) population for both species. The 
electron spectrum is very broad and steadily decreasing with energy. The ion flux is below 
measurable in the 30 to 1000 eV range and then approximately flat to the highest energies measured. 


Day 244.69426 with spacecraft at -5.4 V 



Figure 3. Measured electron and ion spectra at day 244.69426 with chassis at -5.4 V. 
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3.4 Spectra In Sunlight 


Figure 4 shows how the electron spectra can change over time during eclipse and after eclipse 
exit. The lower set of curves were measured during eclipse, when the spacecraft chassis potential 
was -1 to -6 kV. The low energy electrons, which peak at about 10 eV, are secondary electrons 
trapped and accelerated by local electric fields (due to differential charging). When the 
spacecraft exits eclipse, the sunlit surfaces photoemit, generating a large number of low energy 
electrons. Some of these low-energy electrons are also trapped and accelerated by the local 
electric fields. Figure 4 shows enhanced fluxes up to about 60 eV. 

Spectra over 144 minutes on day 272 of 2001 
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Figure 4. Spectra over a 144-minute period on day 272 of 2001, starting in eclipse. The spectra with the higher 
fluxes at lower energies are sunlit. The potential as determined by the ion line ranges between -6339 V and - 
948 V in eclipse and between -718 V and -141 V in sunlight. The line is the flux corresponding to one-count. 

The position and the universal time can be used to determine if the spacecraft is sunlit and the 
incident sun direction. To do this, we implemented the GEOPACK algorithm developed at the 
Goddard Space Flight Center. We confirmed that the algorithm is implemented correctly by 
checking the location of the terminator. All of the spectra checked by hand were correctly 
identified as eclipse or sunlit. 


4 Moments 


For the various datasets, we computed electron and ion density and temperature moments for the 
measured spectra. The results verify that the calculations done at SAIC using the Java software 
using the spin-averaged spectrum give the same results as those done at LANL using the full 
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angular distribution, with the exception of the density of ions with measured energies over 
1 00 eV. In addition to providing verification that the calculations are done correctly, doing these 
calculations allowed us to refine the calculational technique and become familiar with the data. 

The phase space density is derived from the differential flux using the formula 


f = 


2Ee" 


F 


3 6 • • 2 1 

where f is in s m' , E is in eV, and F is in (m s sr eV)" The electron density and temperature 
moments are given by 


47te 2e 


n = 


m 


m 
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where E m is the energy at infinity, E is the geometric mean of the edges of the energy channel, 
and AE is the width of the energy channel. The value of E m ; n for electrons is the maximum of the 
barrier energy (the energy at which the minimum electron count rate between 30 eV and 400 eV 
is observed ) and 30 eV. (There are some additional factors that go into the computation of Emm 
that are described in [Thomsen, et al., 1999].) The value of E m ; n for ions is 124 eV. The 
minimum for electrons is chosen to filter out the secondary and photoelectrons accelerated by 
local electric fields into the MPA. 

For ions, a “low energy” density is also computed in which the sum is over all energy bins below 
E m i n = 124 eV. 

The figures below show the density and temperature moments for electrons and ions in eclipse 
and in sunlight. While only the Sept 01 moments are shown here, the moments using the other 
datasets are similar. 
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LANL computed electron density (cm' 3 ) 


Sept 01 ; Eclipse 



Figure 5. Validity of electron density moment calculational technique for eclipse data. 


Sept 01; Sunlit 



Figure 6. Validity of electron density moment calculational technique for sunlit data. 
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Sept 01 ; Eclipse 



SAIC computed electron temperature (eV) 

Figure 7. Validity of electron temperature moment calculational technique for eclipse data. LANL uses the 
full angular spectrum to compute temperatures parallel and perpendicular to the magnetic field. 


Sept 01; Sunlit 



SAIC computed electron temperature (eV) 


Figure 8. Validity of electron temperature moment calculational technique for sunlit data. LANL uses the full 
angular spectrum to compute temperatures parallel and perpendicular to the magnetic field. 
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Sept 01 ; Eclipse 



Figure 9. Validity of ion density moment calcuiational technique in eclipse. Integral using full angular 
spectrum gives higher densities than spin-averaged spectrum. 


Sept 01 ; Sunlit 



Figure 10. Validity of ion density moment calcuiational technique in sunlight. Integral using full angular 
spectrum gives higher densities than spin-averaged spectrum. 
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Sept 01 ; Eclipse 



SAIC computed ion density (LE) (cm -3 ) 

Figure 11. Validity of low energy ion density moment calculationai technique in eclipse. Spectra for chassis 
potential more negative than -100 V have a low energy ion density of zero. 


Sept 01 ; Sunlit 



SAIC computed ion density (LE) (cm -3 ) 

Figure 12. Validity of low energy ion density moment calculationai technique in sunlight. Spectra for chasiss 
potential more negative than -100 V have a low energy ion density of zero. 
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Sept 01 ; Eclipse 



SAIC computed ion temperature (eV) 


Figure 13. Validity of ion temperature moment calculationai technique in eclipse. Integral using full angular 
spectrum gives lower temperatures than spin-averaged spectrum. 


Sept 01 ; Sunlit 



SAIC computed ion temperature (eV) 


Figure 14. Validity of ion temperature moment calculationai technique in sunlight. Integral using full angular 
spectrum gives lower temperatures than spin-averaged spectrum. 
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The electron moment results computed at SAIC and LANL are in good agreement. This suggests 
that the SAIC implementation of the integrals correctly follows LANL’s approach. There is a 
systematic difference in the ion moments. The density integral using the full angular spectrum 
gives higher densities than the spin- averaged spectrum does. In the temperature moments of 
Figure 13 and Figure 14, the density used in the SAIC computation is the SAIC computed 
density. If the LANL density computed using the full angular spectrum is used instead, Figure 15 
and Figure 16 result. When computed in this way, the spin-averaged temperatures are close to 
the parallel temperatures, while the agreement with the perpendicular temperature is only good 
on average. 


Sept 01 ; Eclipse 



SAIC computed ion temperature (eV) * SAIC den/LANL den 

Figure 15. Validity of ion temperature moment caiculational technique in eclipse using density from full 
angular spectrum. 
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Sept 01 ; Sunlit 



SAIC computed ion temperature (eV) * SAIC den/LANL den 


Figure 16. Validity of ion temperature moment calculational technique in sunlight using density from full 
angular spectrum 


The following figures show relationships between the various moments. For both electrons and 
ions, the temperature varies as the inverse 3/4 power of the density with a great deal of scatter. 
This dependency was determined by plotting the temperature density relation for electrons in 
eclipse and selecting the power that matches the lower curve best. For electrons, there are two 
separate populations that can be present. The higher temperature population is the fresh nightside 
plasma sheet and the lower temperature population is plasma that has drifted around the earth. 
Open triangles are used to distinguish the lower temperature population in several of the figures 
that follow. Points plotted with an open triangle are those that satisfy the relation 

On 0 ' 75 < 300 eVcm -2 ' 25 . 

The ion and electron densities are correlated and comparable. This is consistent with the fact that 
we expect the plasma to be neutral on a length scale comparable with the Debye length 
(hundreds of meters). The ion and electron temperatures are independent. Ion temperatures are 
generally higher and vary over a smaller range. 


17 


Sept 01 ; Eclipse 



Figure 17. Temperature-density relationship for electrons during eclipse. Points plotted with open triangles 
have a value of 9n^' 75 < 300 eVcm “ . 


Sept 01 ; Sunlit 



Figure 18. Temperature-density relation for electrons in sunlit data. Points plotted with open triangles have a 
value of 0n^' 75 < 300 eVcm “ . 
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Sept 01 ; Eclipse 



Figure 19. Temperature-density relationship for ions in eclipse. 


Sept 01 ; Sunlit 



Figure 20. Temperature-density relationship for ions in sunlight. 
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Sept 01 ; Eclipse 



Figure 21. Ion and electron densities are correlated and comparable in eclipse. Points plotted with open 

0 75 —2 25 

triangles have a value of 0n < 300 eVcm for electrons. 


Sept 01 ; Sunlit 



Figure 22. Ion and electron densities are correlated and comparable in sunlight. Points plotted with open 

0 75 —2 25 

triangles have a value of On < 300 eVcm for electrons. 
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Sept 01 ; Eclipse 
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Figure 23. Ion and electron temperatures are independent in eclipse. Ion temperatures are generally higher 

0 75 —2 25 

and vary less. Points plotted with open triangles have a value of 0n < 300 eVcm for electrons. 


Sept 01 ; Sunlit 
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Figure 24. Ion and electron temperatures are independent in sunlight. Ion temperatures are higher and vary 

0 75 —2 25 

less. Points plotted with open triangles have a value of On < 300 eVcm for electrons. 
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LANL has previously reported that there is a relationship between the electron temperature 
calculated as described above and the spacecraft chassis potential. Separate power law fits were 
developed for eclipse and sunlit conditions. For eclipse 

(j) = 1.74xl0“ 5 9 2 - 26 
and for sunlight 

4> = 2.4xlO -7 0 2 ' 51 . 

The following figures show this relationship for the September 01 and March 01 datasets being 
used here. The previous fit, shown by the lines in the figures, is consistent with these datasets. 
The potentials of the low energy electron population points, shown with open triangles, are under 
30 eV and independent of the temperature. 


Sept 01 ; Eclipse 



SAIC computed electron temperature (eV) 


Figure 25. Correlation between chassis potential and electron temperature in eclipse. Points plotted with open 

0 75 —2 25 

triangles have a value of On < 300 eVcm for electrons. 


22 


Sept 96; Sunlit 
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Figure 26. Correlation between chassis potential and electron temperature in sunlight. Points plotted with 

0 75 —2 25 

open triangles have a value of 0n <300eVcm ' for electrons. 


Sept 96; Eclipse 



SAIC computed electron temperature (eV) 


Figure 27. Correlation between chassis potential and electron temperature in eclipse. Points plotted with open 

0 75 —2 25 

triangles have a value of 0n < 300 eVcm for electrons. 
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Sept 96; Sunlit 
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Figure 28. Correlation between chassis potential and electron temperature in sunlight. Points plotted with 

0 75 —2 25 

open triangles have a value of 0n <300eVcm for electrons. 


5 Fluxes 

Figure 29 and Figure 30 show the net electron flux to the spacecraft and its components, incident 
electron flux, secondary electron flux, and backscattered electron flux. The net flux is defined as 

oo 

Flux^ et =-jF e (E)[l-yY e (E)-B(E)]dE= ^ F e (E)(-l + yY e (E) + B(E)) AE 


where the Ys and B are the secondary and backscatter yield functions and y is the fraction of the 
low energy electrons that escape from the spacecraft. The secondary and backscatter yield values 
used are the best available values for carbon graphite and all the low-energy electrons are 
assumed to escape from the spacecraft (y=l). 

The flux is plotted along the x-axis and the potential along the y-axis to show how the potential 
relates to the flux. In order to minimize contamination by secondary and photoelectrons, the 
minimum potential in the sum is given by the same E min as used in the calculations of the 
moments. 

At the higher potentials, the incident electron current goes down as the spacecraft potential goes 
up as a larger fraction of the incident flux is attenuated. 
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Figure 3 1 and Figure 32 show the net ion flux to the spacecraft and its components. Different x- 
axis scales are used for the ion and electron plots. The ion sum is over all energies for which the 
ion flux exceeds twice the one-count flux. 

oo 

Fluxj let = | F i (E)[l + YY i (E)]dE = ^F i (E)(l + yY i (E))AE 
-<t> E 

Based on Figure 29 through Figure 32, it would be difficult to make a case that the potential is a 
function of the incident electron flux, the net electron flux, the incident ion flux, or the net ion 
flux. 


Sept 01 ; Eclipse 



Electron Fluxes (1/(cm 2 s sr)) 


Figure 29. Electron fluxes at all potentials in eclipse. 
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Figure 30. Electron fluxes at all potentials in sunlight. 
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Figure 31. Ion fluxes at all potentials in eclipse. 
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Sept 01; Sunlit 
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Figure 32. Ion fluxes at all potentials in sunlight. 


The net flux to the spacecraft is the sum of the net electron flux and the net ion flux. At 
equilibrium, the total flux to the spacecraft, which consists of the net flux as defined here plus 
any photoemission, is zero. If we have a perfect measurement of the spectrum, perfect 
knowledge of the material properties (secondary and backscatter yields), and there are no 
geometric effects, the computed net flux in eclipse would be zero at all potentials. Figure 33 
shows a small positive net flux with a large amount of scatter for the entire range of potentials. 
On the average, the net flux is 10% of the incident electron flux and greater than the incident ion 
flux, so we need to improve the yield functions and the calculational technique. 

Figure 34 shows the net flux in sunlight. Ignoring all the other contributions to the error in our 
computation of the total flux, the net flux shown is equal in magnitude to the photoemission flux. 
A more energetic environment with higher fluxes charges the spacecraft more. The incident 
electron flux is not attenuated significantly at the lower potentials seen in sunlight charging. 
Therefore the photoemitted flux that is able to escape increases with potential. 

Figure 35 and Figure 36 show the same results for the September 1996 dataset. 
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Sept 01 ; Eclipse 



Fluxes (1/(cm 2 s sr)) 


Figure 33. Net flux in eclipse. A complete spectrum, correct yield functions, and proper accounting for 
suppression of secondaries by barriers would give a value of zero for all potentials. 

Sept 01 ; Sunlit 



Fluxes (1/(cm 2 s sr)) 


Figure 34. Net flux in sunlight. The net flux, as defined here, does NOT include photoemission. The 
photoemission needed for current balance is increasing with potential. 
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Figure 35. Net flux in eclipse. 
Sept 96; Sunlit 



Fluxes (1/(cm 2 s sr)) 


Figure 36. Net flux in sunlight. 


6 Chassis Potential Computation 

We added to our software the capability to compute the potential that we expect given incident 
electron and ion spectra. Comparing the measured and computed potentials provides a measure 
of the quality of the flux computations. 

Conceptually, the computation consists of the following three steps: 

1 . Given the measured potential, compute the ion and electron spectra at infinity. 

2. Using the computed ion and electron spectra at infinity, determine the net current 
(incident, secondaries, backscattered, and photo) to the spacecraft as a function of the 
spacecraft chassis potential. 

3. Search for a unique chassis potential at which the net current is zero. This is the 
computed chassis potential. 

If a computed potential can be found and if the measured and computed potentials are the same 
within the experimental error bar, we conclude that the flux spectra are adequate to compute the 
chassis potential. 

Three items require addressing: 

1. Technique to adjust incident spectra for spacecraft chassis potential. (Required for steps 1 
and 2.) 

2. Calculation of net current as a function of proposed chassis potential. (Step 2) 

3. Technique to search for a unique zero in the net current as function of chassis potential. 
Alternatives to a unique zero. (Step 3) 

4. Size of error bars. (Evaluating results) 

The technique used to adjust the spectra to account for the spacecraft potential is addressed in 
Section 3.2. 

6.1 Net Current 

The net current to the spacecraft in terms of the proposed potential 0' and the measured potential 
0m is given by 
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where F net is the signed average net flux to the spacecraft, (|> m is the measured potential, ())' is the 
proposed potential, the Ys and B are the secondary and backscatter yield functions, y is the 
fraction of the low energy electrons that escape from the spacecraft, f sun iit is the fraction of the 
spacecraft sunlit (about 1/2), F e and Fi are the measured fluxes, and F pho t 0 is the average 
photoemitted flux. 

The electron integral is only strictly correct if the proposed potential is more negative than the 
measured potential, <t>" < (j) m . The portion of the electron spectrum that is not measured, as it is 
does not have enough energy to reach the spacecraft at the measured potential, ideally should be 
included in the integral for the proposed potential. For now, we ignore this discrepancy. 

When using measured fluxes, the integral becomes a sum with the value of E m taken to be the 
energy at the geometric center of the energy bin. Contributions to the ion integral in which the 
measured flux is less than twice the one-count flux are ignored. Contributions to the electron 
integral in which E m < E min are discarded as the flux measurement is dominated by secondary 
and photoelectrons. For both the ion and electron integrals, the lowest energy bin included is the 
one for which the E m value exceeds the lower limit of the integral. Half of the time this is an 
overestimate and half of the time it is an underestimate. 

6.2 Zero In Net Current 

A zero in the net current as a function of potential is searched for between the limits of-1 and 
-10000 V. The current-voltage curve is assumed to have a negative slope over the entire range of 
potential. This is not always true. As a result, the algorithm sometimes fails. There are also cases 
in which the computed net potential does not have a zero in the range of -1 V to -10000 V. The 
algorithm also fails for these cases. When the algorithm fails it returns a value of-1 V. At each 
step, a guess at the potential is made and the net current is computed. A check is made for a net 
current of nearly zero (absolute value under 10’ 6 A m" 2 ) or a change in sign of the net current 
between successive guesses. When a zero in the curve is found, a linearly interpolated value is 
returned. The initial guess is the measured potential and each successive guess is the previous 
guess times or divided by 1.01 depending on the sign of the net current. Presently, the zero is 
assumed to be unique and the possibility of multiple zeros remains to be explored. 

Often the above algorithm fails to provide a solution. We also use another procedure in which we 
assume that the net current is always positive and look for the potential with the minimum value 
of net current. This procedure is only applied when the first fails. 
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6.3 Error Bars 


The measured potential has a sizeable error bar, which also leads to error bars in the net current 
calculations that depend on the measured potential. 

Throughout most of the energy range of the MPA instrument, the energy of the upper edge of the 
energy bin is approximately 1.31 times the energy of the lower edge. Neighboring energy bins can 
also respond to an incident particle. The measured potential is the geometric center of the lowest 
energy bin with a significant number of counts. (The actual algorithm is much more complex than 
this, but for present purposes this simplification is adequate.) Therefore, the actual potential is only 
measured within a factor of 1 . 1 5 . c)) m /1 . 1 5 < 4> actual <1.1 5(j) m . Additional errors are introduced in the 
computation of the net current using the measured potential, rather than the actual potential. And 
finally, we are treating the spacecraft as a uniform sphere, introducing significant errors. Somewhat 
arbitrarily, we take agreement within a factor of 1 .5 to be good agreement. 

6.4 Results 

Figure 37 and Figure 38 compare the measured and computed potential for the September 2001 
eclipse dataset. In Figure 37, we look for a zero in the net flux. Loosening the criteria to search 
for a minimum in the flux rather than a zero, increases the number of cases in which a solution 
can be found as most of the solutions are good solutions. The poor potential estimate shown in 
these figures tells us that we need to improve our technique of computing the net current. 

Sept 01 ; Eclipse 



Sept 01 ; Eclipse 



Number no solution 


Figure 37. For Sept 01 eclipse data, measured potential as a function of the potential computed from the zero 
in the net flux. The lines are y = 1.5 x and y = x/1.5. The bar chart at the right shows the number of spectra 
for which no solution was found for each value of measured potential. These points are represented in the plot 
on the left by a single diamond along the left hand axis. 
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Figure 38. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux. The lines are y = 1.5 x and y = x/1.5. The bar chart at the right shows the number of 
spectra for which no solution was found for each value of measured potential. 


7 Calculational Techniques 

7 . 1 Yield Functions 

The computations of fluxes shown above require knowledge of the average yield properties of 
the spacecraft surfaces. We are treating the spacecraft as if it was covered entirely by a single 
average material. The properties we choose to use for the first set of calculations described here 
are for graphite and are shown in Table 1 . The electron-generated secondary yield properties are 
chosen to fit the curve in Barnett, et al. [1977], which tabulates data from Holzland and Jacobi, 
[1969] and Bruining [1938]. The ion-generated secondary yield properties are chosen to be 
consistent with the curve for incident Hydrogen ions in Barnett, et al. [1977], which tabulates 
data from Large and Whitlock, [1962]. The effective atomic number, which is used to compute 
the backscatter, is the default value for graphite in the SEE Spacecraft Charging Handbook. 
Figures 39-42 show the yield properties of this material, labeled “Carbon,” and a selection of 
materials from the default materials supplied with the SEE Spacecraft Charging Handbook. We 
also reviewed preliminary results from measurements of yield properties made at Utah State. The 
yield curves derived from the default properties of the Gold, Kapton, OSR, and Solar Cell 
materials included in the SEE Spacecraft Charging Handbook are representative of the range of 
spacecraft materials. Flux and potential calculations using the properties of these four materials 
are shown in Figures 43-50. All of the other sets of material properties give much larger positive 
net current than “Carbon” does in the hundreds through thousands of volts range of spacecraft 
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potential. The estimated potentials using the OSR and Solar cell material properties are 
consistently higher than the measured potentials. Using Gold material properties tends to 
overestimate the potential, and using Kapton material properties sometimes overestimates it and 
sometimes underestimates it. The low secondary yield peak at low incident electron energy 
represents the actual average spacecraft material the best. The rest of the calculations shown in 
this report all use the properties shown in Table 1. 


Table 1. Material properties used average material “Carbon” 


Property 

Value 

Backscatter 

Effective atomic no. 

4.5 

Secondary electron yield due to incident electrons 

Max sec. yield 

0.93 

Energy of peak 

0.28 keV 

Range Param 1 

180 

Range Exp. 1 

0.45 

Range Param 2 

312 

Range Exp. 2 

1.75 

Secondary electron yield due to incident ions 

Sec. Yield of 1 keV proton 

0.455 

Proton energy 

80 keV 
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Figure 39 Secondary electron yield from incident electron impact for materials defined in the SEE Spacecraft 
Charging Handbook and the Carbon material used here. 



Incident electron energy (eV) 


Figure 40. Backscattered electron yield from incident electron impact for materials defined in the SEE 
Spacecraft Charging Handbook and the Carbon material used here. 
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Figure 41 Sum of secondary and backscattered electron yields from incident electron impact for materials 
defined in the SEE Spacecraft Charging Handbook and the Carbon material used here. 



Incident ion energy (eV) 


Figure 42 Secondary electron yields from incident ion impact for materials defined in the SEE Spacecraft 
Charging Handbook and the Carbon material used here. 
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Sept 01; Eclipse; Gold spacecraft 



Fluxes (1/(cm 2 s sr)) 


Figure 43. Net flux in eclipse using material properties of Gold as defined in the SEE Spacecraft Charging 
Handbook. 


Sept 01; Eclipse; Gold spacecraft 



Figure 44. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux using material properties of Gold as defined in the SEE Spacecraft Charging 
Handbook. 
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Sept 01; Eclipse; Kapton spacecraft 



Fluxes (1/(cm 2 s sr)) 


Figure 45. Net flux in eclipse using material properties of Kapton as defined in the SEE Spacecraft Charging 
Handbook. 


Sept 01; Eclipse; Kapton spacecraft 



Figure 46. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux using material properties of Kapton as defined in the SEE Spacecraft Charging 
Handbook. 
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Sept 01; Eclipse; Solar cell spacecraft 
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Figure 47. Net flux in eclipse using material properties of Solar cell coverglass material as defined in the SEE 
Spacecraft Charging Handbook. 


Sept 01 ; Eclipse; Solar cell spacecraft 



Figure 48. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux using material properties of Solar cell coverglass material as defined in the SEE 
Spacecraft Charging Handbook 
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Sept 01; Eclipse; OSR spacecraft 
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Figure 49. Net flux in eclipse using material properties of OSR as defined in the SEE Spacecraft Charging 
Handbook. 


Sept 01; Eclipse; OSR spacecraft 



Figure 50. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux using material properties of OSR as defined in the SEE Spacecraft Charging 
Handbook. 
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7.2 High Energy Extrapolation 


The MPA measurements only extend up to 47 keV for electrons and 49 keV for ions. We were 
interested in exploring the possibility that we might be able to obtain better results by extending 
the integral beyond the measurements. We do this by fitting the last five points of the spectrum 
to a Kappa distribution function. We then add additional points to the sum up to an energy of 10 5 
eV ± cf), where the “+” is for the electron spectrum and for the ion spectrum. Each point is for 
an energy of 1 .3 1 times the previous point. No checking for quality of fit is done. Figure 5 1 and 
Figure 52 show the resulting fluxes. This approach does not appear to be useful and has not been 
pursued further. 


Sept 01; Eclipse; Extended integral 



Fluxes (1/(cm 2 s sr)) 


Figure 51. Net flux in eclipse including contribution from extrapolated fluxes. 
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Sept 01; Sunlit; Extended integrals 



Fluxes (1/(cm 2 s sr)) 


Figure 52. Net flux in sunlight, including contribution from extrapolated fluxes. 


7.3 Low Energy Electron Cutoff 

Another aspect of the integral that was explored is the lower energy cutoff in the electron 
spectrum. When the spacecraft is charged, the low energy part of the electron spectrum is usually 
dominated by low energy secondary and photo electrons trapped by differential potentials near 
the detector. The lower limit of the electron integrals needs to be as low as possible while high 
enough to eliminate the trapped electrons. LANL computes this energy by searching for the 
minimum in the count rate in the energy bins between 30 eV and 400 eV. [Thomsen, et al, 1999] 
The energy at the center of this bin is reported as the barrier. Figure 53 shows this barrier as a 
function of potential for the September 2001 dataset sunlit points. For eclipse points, the “LANL 
barrier” is always 34. 18 eV. For points with a barrier of 34. 18, we propose an alternative barrier 
of half the potential. These points are shown as triangles in Figure 53. Figure 54 shows the fluxes 
in eclipse using the alternative barrier rather than the LANL barrier. The only significant 
differences are in the low potential points. When the potential is computed, Figure 55, many of 
the low potential points falls on within the 1.5 range of a perfect fit. Figure 56 shows the sunlit 
fluxes. We also tried some other approaches to determining the appropriate value to use for the 
cutoff. Nothing we tried works better than these two values. The rest of the calculations in this 
report use the alternative barrier. 
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Sept 01 ; Sunlit 
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Figure 53. Minimum energy used in the calculations up to this point (LANL barrier) and an alternative as a 
function of potential for sunlit calculations. 


Sept 01; Eclipse; Alternative minimum e energy 
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Figure 54. Net flux in eclipse using “alternative barrier” as the minimum electron energy. 
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Sept 01; Eclipse; Alternative minimum e energy 
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Figure 55. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux using “alternative barrier” as the minimum electron energy. 


Sept 01; Sunlit; Alternative minimum e energy 
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Figure 56. Net flux in sunlight using “alternative barrier” as the minimum electron energy. 
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7.4 Fraction Of Low Energy Electrons Escaping 


Another factor that is not included in the above flux calculations is that not all secondary and 
photo electrons escape the spacecraft. (Backscattered electrons have higher energies and are 
unlikely to be trapped.) Let y be the fraction that escape to infinity. Then we have that the total 
current, which on the timescale of MPA measurements is always equal to zero, is given by 

Total flux — 0 — — F e + Fj + + yF sec + yfphoto • 

For each energy, we then have an expression for y in terms of quantities that we have been 
calculating. 


y _ Fc fimek 

' p p 

r sec ~ r r photo 

The low energy flux to the MPA, which we have been carefully keeping out of our flux integrals, 
consists of low energy electrons that do not escape. The fraction of the created low energy flux 
that is measured by the MPA is given by 

o _ measured low energy flux 

^ p p 

1 see ~ r r photo 

If the MPA is located at an “average” point on the spacecraft, y + (3 ~ 1 . These quantities are 
shown for eclipse periods and sunlit periods in Figure 57 and Figure 58 respectively. In eclipse, 
the average value of y + (3 computed in this way is 0.92, slightly under 1.0. In eclipse, the 
measured low energy flux is a useful representation of the average low-energy return flux. The 
value of y is about 0.8. As expected, this approach does not work as well in sunlight. There is a 
great deal of scatter in the points. A photoemission flux of 2.9 x 10 (cm“ s sr)" (2.9 x 10' A m" 
2 ) is needed to make the average value of y + (3 equal 1.0. This is approximately half of the value 
given in Feuerbacher and Fitton [1972] for photoemission flux from graphite in sunlight. With 
this value for photoemission, the average value of y is 0.28 and the average value of (3 is 0.72. 

If we assume that only 81% of the secondary electrons escape, the net fluxes look as shown in 
Figure 59. The net flux is essentially zero over the entire range. The estimated potential is shown 
in Figure 60. 65% of the points have a solution within a factor of 1.5, 87% of the points have a 
solution within a factor of 3, and 8% cannot find a solution. Almost all of the points for which no 
solution can be found have a potential below 30 V. 
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Sept 01; Eclipse; Alternative minimum e energy 



Figure 57. Measures of escaping and returning low energy fluxes in eclipse. 


Sept 01; Sunlit; Alternative minimum e energy 



Figure 58. Measures of escaping and returning low energy fluxes in sunlight. 
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Sept 01; Eclipse; 81% low energy electrons escape 



Fluxes (1/(cm 2 s sr)) 


Figure 59. Net flux in eclipse with alternative minimum electron energy and 81% escape fraction. 


Sept 01; Eclipse; 81% low energy electrons escape 
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Figure 60. For Sept 01, eclipse data, measured potential as a function of the potential computed from the 
minimum in the net flux with alternative minimum electron energy and 81% escape fraction. 
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7.5 Modeling Sunlit Charging 


In sunlight, photoemission that escapes from the spacecraft contributes to the net current. Above, 
we wrote the photoemission contribution to the net flux as y p ((|/) f sunlit F photo (([)') . In order to 

estimate this term, we first look at the net flux, Figure 61, and estimated potential, Figure 62, 
without this term. The net flux is generally negative and is increasing negative as the potential 
increases. The way to understand this dependence is as follows. The amount of photoemission 
created is independent of the potential, but the barrier that develops with time depends on the 
differential potential. A more energetic environment with higher fluxes charges the spacecraft 
more. The incident electron flux is not attenuated significantly at the under kilovolt potentials 
seen in the sunlight charging. Therefore the photoemitted flux that is able to escape increases 
with potential. 

The estimated potential ignoring photoemission is an overestimate as the computed net flux is 
negative, rather than near zero as it would be if we accounted for photoemission. 

In order to estimate the photoemission term, we plot the average net flux as a function of 
potential, Figure 63. We fit a straight line to the portion of the semi-log plot that looks linear. 

The fit gives J photo (([)) = 8xl0 7 (0.171n(|c|)|)-0.38) in units of l/(cm 2 s sr). 8 x 10 7 /(cm 2 s sr) is 

equivalent to 8 x 10' Amps/m . Feuerbacher and Fitton [1972] give a value of 7.2 x 10' 
Amps/m 2 for the average photoemission rate of graphite in sunlight. Assuming a fraction sunlit 
of about 1/2, the coefficient is about 5 times smaller than the photoemission flux from graphite. 
We estimate the photoemission term to be given by 


Yp(0) f suniit F photo(^) = 8xlo7xmin ( 1 ’ max (°’( O - 171n (H) _O - 38 ))) 


At potentials less negative than -9.4 V, the escaping photoemitted current is zero. At potentials 
more negative than -3.35 kV, the escaping photoemitted current is 8 x 10 /(cm s sr), about 20% 
of the available photoemitted current. 

Figure 65 shows the resulting potentials as a function of the estimated potential. 48% of the 
points are within 1 .5 of the measured potential and 81% of the points are within a factor of 3 of 
the measured potential. This includes 21% for which the algorithm correctly identified potentials 
near or below 1 V. Overall, this procedure provides a moderately quantitative prediction of 
chassis potential in sunlight. 

That the results are less well behaved than the eclipse data is not surprising. The development of 
potential barriers and the resulting amount of photocurrent that escapes is, in general, strongly 
dependent on the spacecraft shape and orientation with respect to the sun and the time in which 
the spacecraft has been charging along with the incident electron and ion spectra. 
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Sept 01; Sunlit; 81% low energy electrons escape 
Photoemission not included 



Fluxes (1/(cm 2 s sr)) 


Figure 61. Net flux in sunlight with alternative minimum electron energy, 81% escape fraction, and ignoring 
photoemission. 


Sept 01; Sunlit; 81% secondary electrons escape 
Photoemission not included 



Figure 62. For Sept 01, sunlit data, assuming 81% of the secondary electrons escape, and ignoring 
photoemission, measured potential as a function of potential computed from the minimum in the net flux. The 
lines are y = 1.5 x and y = x/1.5. 


49 


Sept 01; Sunlit; 81% of secondaries escape 
Not including photoemission 

Measured potential (V) 


1 10 100 1000 10000 



Figure 63. Average net flux (ignoring photoemission) as a function of potential for the Sept 01, sunlit dataset. 
The line is a fit to the average net flux for potentials between -30 and -948 V. 


Sept 01 ; Sunlit 
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Figure 64. Net flux in sunlight with alternative minimum electron energy, 81% escape fraction, and including 
photoemission. 
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Sept 01; Sunlit; 81% secondary electrons escape; 
include photoemission 



Figure 65. For Sept 01, sunlit data, assuming 81% of the secondary electrons escape, and including photoemission, 
measured potential as a function of potential computed from the minimum in the net flux. The lines are y = 1.5 x and 
y = x/1.5. The bar chart at the right shows the number of spectra for which no solution was found for each value of 
measured potential. These points are represented in the plot on the left by a single diamond along the left hand axis. 

7.6 Alternative Predictors 

In evaluating our ability to predict potential from the measured spectra, we need to compare our 
approach with other approaches. 

LANL has previously reported that there is a relationship between the electron temperature 
calculated as described in Section 0 and the spacecraft chassis potential. Separate power law fits 
were developed for eclipse and sunlit conditions. For eclipse 

4> = 1.74xlO -5 0 2 ' 26 
and for sunlight 

(j) = 2.4xlO“ 7 0 2 ' 51 . 

The relationship only applies if the flux is above a threshold value. In eclipse, this relationship, 
shown in Figure 55, is as good a predictor as ours above a kilovolt, but gives consistently too 
high a potential below that. It does not work when the measured potential is below about 30 V. It 
may predict potentials of the order of 100 V. 

In sunlight, this relationship, shown in Figure 67, doesn’t work quite as well as the one we 
developed, shown in Figure 65. 30% of the points are within a factor of 1.5 and 85% are within a 
factor of 3. It underestimates the potential by more than a factor of 1 .5 for more than half of the 
spectra. It also often predicts a significant potential when the measured potential is under 1 V (a 
large peak in the lowest energy ion channel). 
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Sept 01 ; Eclipse 



Figure 66. For Sept 01, eclipse data, measured potential as a function of (|) = 1 .74x10 5 0 22>6 for eclipse Sept 
01 dataset. 


Sept 01; Sunlit 



Figure 67. For Sept 01, sunlit data, measured potential as a function of <|) = 2.4x10 0 2 M for sunlit Sept 01 
dataset. 
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We also examined the incident electron flux in the higher energy channels as possible predictors 
of the potential. We examined single channels and all channels above a specific energy. The best 
predictor for the September 2001, eclipse dataset is the sum of the bins from 9123 V and above. 
The best fit to this dataset is 

^estimated = 1-^59 xl0“ 12 ^ (measured flux x bin width) - "' . 

E>9123 

This approach always provides a potential estimate. It works best at the highest potentials. It is 
not as good a predictor as the temperature. 


Sept 01 ; Eclipse 



Figure 68. For Sept 01, eclipse data, measured potential as a function of 
1.359xl0~ 12 ^ (measured flux x bin width)" 03 for sunlit Sept 01 dataset. 

E>9123 


8 Comparison Of Solar Maximum And Solar Minimum 

As shown in Figures 69-72, the fluxes and resulting charging during charging periods are the 
same at solar minimum as at solar maximum. What varies is the frequency of high charging 
environments, with less frequent charging events at solar minimum. 
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Sept 96; Eclipse; 81% low energy electrons escape 



Figure 69. For Sept 96, net flux in eclipse with alternative minimum electron energy and 81% escape fraction. 
(Compare with Figure 59.) 


Sept 96; Eclipse; 81% low energy electrons escape 



Sept 01 ; Eclipse 



Figure 70. For Sept 96, eclipse data, assuming 81% of the secondary electrons escape, measured potential as a 
function of potential computed from the minimum in the net flux. (Compare with Figure 60.) 
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Sept 96; Sunlit; 81% low energy electrons escape 



Figure 71. For Sept 96, net flux in sunlight with alternative minimum electron energy, 81% escape fraction, 
and including photoemission. (Compare with Figure 64.) 


Sept 96; Sunlit; 81% low energy electrons escape 



Potential from zero in flux (V) 


Figure 72. For Sept 96, sunlit data, assuming 81% of the secondary electrons escape, and including 
photoemission, measured potential as a function of potential computed from the minimum in the net flux. 
(Compare with Figure 65.) 
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9 Comparison With Data From Another LANL Spacecraft 


To further check the validity of our calculations, we did a single calculation using a dataset from 
another of the LANL spacecraft. The results are shown in Figure 73 and Figure 74. 

LANL 97; Sept 01 ; Eclipse; 81 % low energy electrons escape 



Fluxes (1/(cm 2 s sr)) 


Figure 73. For LANL 97 spacecraft, net flux in eclipse with alternative minimum electron energy and 81% 
escape fraction. (Compare with Figure 59.) 
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Figure 74. For Sept 01, eclipse data from LANL 97 spacecraft, assuming 81% of the secondary electrons 
escape, measured potential as a function of potential computed from the minimum in the net flux. (Compare 
with Figure 60.) 
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10 Fitting Techniques And Results 


The goal of this research is to better understand the plasma environment that causes spacecraft 
charging and determine a relatively simple characterization of the full particle distributions that 
yields an accurate prediction of the observed charging under a wide variety of conditions. The 
computer codes that compute spacecraft charging in a tenuous plasma environment 
(NASCAP/GEO [ Katz, et al, 1979], SEE Spacecraft Charging Handbook [Katz, et al, 2000], 
and NASCAP-2K [Davis, et al, 2002]) use a Maxwellian distribution function. Originally the 
Maxwellian distribution function was chosen for its simplicity and no one has ever established its 
adequacy or inadequacy. In the late 70s and early 80s, when NASCAP/GEO was written, some 
calculations were attempted with measured spectra. The high noise level and low resolution of 
the measured spectra made the solutions unreliable, so Maxwellian fits are used today instead. 
We added to our flux computation code the ability to compute some analytic functional fits to the 
spectra. We can then use fit spectra in our calculations of moments of the distribution, the fluxes, 
and the computed potential. 


The proposed functional forms are appropriate for the flux at infinity. Before fitting the 
measured fluxes to a functional form, the flux is multiplied by the orbit limited l-cf)/E factor and 
the energy is shifted by the potential. 


F-(E.) = F-(E„±^) = F„(Ej 
Six different functional forms have been implemented. 
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F lux Quadratic = ^2 + BE + C 

Flux Exponential (E) = A exp (BE) 
Flux PowerLaw (E) = AE N 


The Kappa distribution was selected as Christon, et al [1989] have shown that it provides a good 
fit to the quiescent plasma sheet environment at greater than 12 Re (Geosynchronous is at 6.6 
R e .). The Kappa function has the shape of a Maxwellian at low energies and a power law at high 
energies, providing a high-energy tail to the distribution. The active magnetosphere is more 
complex [Christon, et al, 1991], We are looking for a simple model of the environment so we use 
a single Kappa distribution. 


The Maxwellian, Double Maxwellian, and Kappa distribution functions were written for the 
purpose of fitting the flux. The exponential and power law fitting procedures were written to 
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assist in the extrapolation of the measured data to higher energies. The quadratic fitting 
procedure was developed for finding the minimum in one of the procedures we tried (and 
discarded) to determine the best value to use as the low energy cutoff in the electron integrals. 

Two least-square fitting procedures were developed. The power law, exponential, and 
Maxwellian functions are fit by taking logarithms of the flux and the energy, as appropriate, and 
computing the best-fit straight line. Each point is weighted by the bin width. The fitting 
procedure (lfit from Numerical Recipes [Press, et al. 1992]) finds values for a and b that 
minimize the expression 

Z( A i(yi _a_bx i )) 2 

i 

where 


Ai 

bin width, 

y; 

logarithm of the flux at infinity (power law and exponential) or the difference 
between the logarithm of the flux at infinity and the logarithm of the energy 
(Maxwellian), 

Xi 

energy (exponential and Maxwellian) or the logarithm of the energy (power law). 


The Double Maxwellian, Kappa, and Quadratic functions are fit using the Levenberg-Marquardt 
method for nonlinear least-squares fitting. The algorithm and an implementation in C are 
described in Numerical Recipes [Press, et al. 1992], Again the points are weighted by the bin 
width. The fitting procedure finds the vector a that minimizes the expression 

XfA^y.-ffa.x ,))) 2 

i 

where 


Ai 

bin width, 

Yi 

flux at infinity, 

Xi 

energy 

f 

functional form 

a 

vector of the constants in the functional form. For a quadratic, a = (A, B, C). For a 
Kappa distribution, a = (A, E 0 , k). For a double Maxwellian distribution, a = (ni, 0i, 
n,*n 2 , 0i*e 2 ). 


Both approaches are used to compute a single Maxwellian fit. The double Maxwellian fitting 
procedure can ignore the third and fourth parameters. We refer to a fit to a single Maxwellian 
using the first approach as a Maxwellian fit and a fit to a single Maxwellian using the second 
approach as a Single Maxwellian fit. 

An alternative fitting procedure in which moments of the distribution are strictly preserved was 
proposed but not implemented. 
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In fitting ion spectra, the portion of the spectrum fit is from the highest energy bin through the 
energy bin just above the potential excluding energy bins in which the flux is below twice the one- 
count equivalent flux. This avoids all energy bins with possibly misleading count rates. 

In fitting electron spectra, the portion of the spectrum fit is the range over which the flux integrals 
are done. This avoids all energy bins with possibly misleading count rates. 

Figures 75 and 76 show fits of a charging spectrum and non-charging spectrum to various 
functional forms. 

Moments and fluxes can be computed from fit distributions in exactly the same way as they are 
computed from the measured distributions. Results are shown in Figures 77-108. With fit 
distributions, it is possible to include energy bins that are not included in the integrals over the 
measured distribution, such as ion fluxes below twice the one-count flux equivalent, electron 
energy bins below the estimated potential barrier, and energy bins above the maximum of the 
instrument, 47 keV. We have the ability to use either approach. Figures 77-92 use the entire energy 
range and Figures 93-108 use only those energy bins used in the integrals over the measured 
distributions. The maximum energy chosen can affect the results a great deal, particularly when 
fitting the ions with a Kappa distribution, as the knee of the ion spectrum is often above 47 keV. 
Presently, the upper limit of the net flux integrals is about 100 keV. 

Figures 109-124 show electron and ion component fluxes, net fluxes, and computed potentials with 
the fluxes computed from fits. These figures should be compared with Figures 29, 31, 59 and 60 
computed using the measured fluxes directly. In these figures, the electrons and ions use the same 
analytic fomi for the fit distribution function. The figures using the Maxwellian, Single 
Maxwellian, and Double Maxwellian fits are similar with no improvement in predictive ability 
using the additional two parameters of the Double Maxwellian. The Kappa distribution provides a 
very poor estimate of the ion fluxes. 

Figures 125-136 show the measured and computed potentials using the various combinations of 
distribution function analytic forms. The accuracy is summarized in Table 2. Using Maxwellian 
fits is never as accurate as using the measured fluxes directly. However, using a Kappa distribution 
for the electrons and either single Maxwellian for the ions gives results similar to those obtained 
from the measured fluxes directly. 

In a number of the measured potential versus computed potential figures, there are a handful of 
points in the lower right hand comer. These points have measured potentials on the order of 10 V 
and computed potentials on the order of a few hundred volts. For several of these cases, the 
temperature for the Maxwellian fit to the ion spectrum is of the order of a thousand and there is a 
peak in the measured spectrum around 10 V. Upon careful examination it appear that the fitting 
procedure is missing the low energy peak as there are far more points in the higher energy portion 
of the spectrum. The fits would provide more accurate potential calculations if additional 
intelligence to look for this feature were added to the fitting procedure. 

Table 3 compares the results including all energy bins in the flux integrals with results using only 
those bins used in the sums computed using the measured fluxes directly. Neither approach 
provides more consistently accurate potential computations. 
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Fits to charging spectrum 



Figure 75. A sample charging flux spectrum and various fits to it using different functional forms. 


Fits to non-charging spectrum 



Figure 76. A sample non-charging flux spectrum and various fits to it using different functional forms. 
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Sept 01; Eclipse; Maxwellian fits 



Figure 77. Electron density computed from data compared with electron density computed from fit to 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 


Sept 01; Eclipse; Maxwellian fits 



Electron temperature computed from fit (eV) 

Figure 78. Electron temperature computed from data compared with electron temperature computed from fit 
to Maxwellian. Integral over fit extends from 1 eV to 100 keV. 
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Sept 01; Eclipse; Maxwellian fits 



Ion density computed from fit (cm' 3 ) 


Figure 79. Ion density computed from data compared with ion density computed from fit to Maxwellian. 
Integral over fit extends from 1 eV to 100 keV. 


Sept 01; Eclipse; Maxwellian fits 



Ion temperature computed from fit (eV) 


Figure 80. Ion temperature computed from data compared with ion temperature computed from fit to 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 


62 


Sept 01; Eclipse; Single Maxwellian fits 



Figure 81. Electron density computed from data compared with electron density computed from fit to Single 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 

Sept 01 ; Eclipse; Single Maxwellian fits 



Electron temperature computed from fit (eV) 

Figure 82. Electron temperature computed from data compared with electron temperature computed from fit 
to Single Maxwellian. Integral over fit extends from 1 eV to 100 keV. 
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Sept 01; Eclipse; Single Maxwellian fits 



Ion density computed from fit (cm' 3 ) 


Figure 83. Ion density computed from data compared with ion density computed from fit to Single 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 


Sept 01; Eclipse; Single Maxwellian fits 



Ion temperature computed from fit (eV) 


Figure 84. Ion temperature computed from data compared with ion temperature computed from fit to Single 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 


64 


Sept 01; Eclipse; Double Maxwellian fits 



Figure 85. Electron density computed from data compared with electron density computed from fit to Double 
Maxwellian. Integral over fit extends from I eV to 100 keV. 


Sept 01; Eclipse; Double Maxwellian fits 



Electron temperature computed from fit (eV) 

Figure 86. Electron temperature computed from data compared with electron temperature computed from fit 
to Double Maxwellian. Integral over fit extends from 1 eV to 100 keV. 
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Sept 01; Eclipse; Double Maxwellian fits 



Ion density computed from fit (cm' 3 ) 


Figure 87. Ion density computed from data compared with ion density computed from fit to Double 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 


Sept 01; Eclipse; Double Maxwellian fits 



Ion temperature computed from fit (eV) 


Figure 88. Ion temperature computed from data compared with ion temperature computed from fit to Double 
Maxwellian. Integral over fit extends from 1 eV to 100 keV. 
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Sept 01; Eclipse; Kappa fits 



Figure 89. Electron density computed from data compared with electron density computed from fit to Kappa. 
Integral over fit extends from 1 eV to 100 keV. 


Sept 01; Eclipse; Kappa fits 



Electron temperature computed from fit (eV) 

Figure 90. Electron temperature computed from data compared with electron temperature computed from fit 
to Kappa. Integral over fit extends from 1 eV to 100 keV. 
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Sept 01; Eclipse; Kappa fits 



Figure 91. Ion density computed from data compared with ion density computed from fit to Kappa. Integral 
over fit extends from 1 eV to 100 keV. 


Sept 01; Eclipse; Kappa fits 



Ion temperature computed from fit (eV) 


Figure 92. Ion temperature computed from data compared with ion temperature computed from fit to 
Kappa. Integral over fit extends from 1 eV to 100 keV. 
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Sept 01; Eclipse; Maxwellian fits 



Figure 93. Electron density computed from data compared with electron density computed from fit to 
Maxwellian. Integral over fit only includes energy bins used to create fit. 

Sept 01 ; Eclipse; Maxwellian fits 



Electron temperature computed from fit (eV) 

Figure 94. Electron temperature computed from data compared with electron temperature computed from fit 
to Maxwellian. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Maxwellian fits 



Ion density computed from fit (cm' 3 ) 


Figure 95. Ion density computed from data compared with ion density computed from fit to Maxwellian. 
Integral over fit only includes energy bins used to create fit. 


Sept 01; Eclipse; Maxwellian fits 



Ion temperature computed from fit (eV) 


Figure 96. Ion temperature computed from data compared with ion temperature computed from fit to 
Maxwellian. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Single Maxwellian fits 



Figure 97. Electron density computed from data compared with electron density computed from fit to Single 
Maxwellian. Integral over fit only includes energy bins used to create fit. 

Sept 01 ; Eclipse; Single Maxwellian fits 



Electron temperature computed from fit (eV) 

Figure 98. Electron temperature computed from data compared with electron temperature computed from fit 
to Single Maxwellian. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Single Maxwellian fits 



Ion density computed from fit (cm' 3 ) 


Figure 99. Ion density computed from data compared with ion density computed from fit to Single 
Maxwellian. Integral over fit only includes energy bins used to create fit. 


Sept 01; Eclipse; Single Maxwellian fits 



Ion temperature computed from fit (eV) 


Figure 100. Ion temperature computed from data compared with ion temperature computed from fit to Single 
Maxwellian. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Double Maxwellian fits 



Electron density computed from fit (cm -3 ) 


Figure 101. Electron density computed from data compared with electron density computed from fit to 
Double Maxwellian. Integral over fit only includes energy bins used to create fit. 

Sept 01; Eclipse; Double Maxwellian fits 



Electron temperature computed from fit (eV) 

Figure 102. Electron temperature computed from data compared with electron temperature computed from 
fit to Double Maxwellian. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Double Maxwellian fits 



Ion density computed from fit (cm' 3 ) 


Figure 103. Ion density computed from data compared with ion density computed from fit to Double 
Maxwellian. Integral over fit only includes energy bins used to create fit. 


Sept 01; Eclipse; Double Maxwellian fits 



Ion temperature computed from fit (eV) 


Figure 104. Ion temperature computed from data compared with ion temperature computed from fit to 
Double Maxwellian. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Kappa fits 



Figure 105. Electron density computed from data compared with electron density computed from fit to 
Kappa. Integral over fit only includes energy bins used to create lit. 

Sept 01; Eclipse; Kappa fits 



Electron temperature computed from fit (eV) 

Figure 106. Electron temperature computed from data compared with electron temperature computed from 
fit to Kappa. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Kappa fits 



Ion density computed from fit (cm -3 ) 


Figure 107. Ion density computed from data compared with ion density computed from fit to Kappa. Integral 
over fit only includes energy bins used to create fit. 

Sept 01; Eclipse; Kappa fits 
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Figure 108. Ion temperature computed from data compared with ion temperature computed from fit to 
Kappa. Integral over fit only includes energy bins used to create fit. 
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Sept 01; Eclipse; Maxwellian fits 
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Figure 109. Electron fluxes in eclipse computed from Maxwellian fit at all potentials. Integral extends from 1 
eV to 100 keV. 


Sept 01; Eclipse; Maxwellian fits 
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Figure 110. Ion fluxes in eclipse computed from Maxwellian fit at all potentials. Integral extends from 1 eV to 
100 keV. 
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Sept 01 ; Eclipse; Max & Max fits 



Fluxes (1/(cm 2 s sr)) 


Figure 111. Net flux in eclipse computed from Maxwellian fit at all potentials. Integral extends from 1 eV to 
100 keV. 


Sept 01 ; Eclipse; Max & Max fits 



Figure 112. Measured potential as a function of the potential computed from the minimum in the net flux, 
where the fluxes are computed from Maxwellian fits to the electron and ion fluxes. 
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Sept 01; Eclipse; Single Maxwellian fits 



Electron Fluxes (1/(cm 2 s sr)) 


Figure 113. Electron fluxes in eclipse computed from Single Maxwellian fit at all potentials. Integral extends 
from 1 eV to 100 keV. 


Sept 01; Eclipse; Single Maxwellian fits 
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Figure 114. Ion fluxes in eclipse computed from Single Maxwellian fit at all potentials. Integral extends from 
1 eV to 100 keV. 



□ incident ion 

□ secondary ion 
o net ion 

Reference 


79 


Sept 01 ; Eclipse; Single Max & Single Max fits 
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Figure 115. Net fluxes in eclipse computed from Single Maxwellian fit at all potentials. Integral extends from 
1 eV to 100 keV. 


Sept 01 ; Eclipse; Single Max & Single Max fits 
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Figure 116. Measured potential as a function of the potential computed from the minimum in the net flux, 
where the fluxes are computed from Single Maxwellian fits to the electron and ion fluxes. 
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Sept 01; Eclipse; Double Maxwellian fits 



Electron Fluxes (1/(cm 2 s sr)) 


Figure 117. Electron fluxes in eclipse computed from Double Maxwellian fit at all potentials. Integral extends 
from 1 eV to 100 keV. 


Sept 01; Eclipse; Double Maxwellian fits 



Ion Fluxes (1/(cm 2 s sr)) 


Figure 118. Ion fluxes in eclipse computed from Double Maxwellian fit at all potentials. Integral extends from 
1 eV to 100 keV. 
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Sept 01 ; Eclipse; Double Max & Double Max fits 



Fluxes (1/(cm 2 s sr)) 


Figure 119. Net fluxes in eclipse computed from Double Maxwellian fit at all potentials. Integral extends from 
1 eV to 100 keV. 
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Figure 120. Measured potential as a function of the potential computed from the minimum in the net flux, 
where the fluxes are computed from Double Maxwellian fits to the electron and ion fluxes. 
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Sept 01; Eclipse; Kappa fits 
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Figure 121. Electron fluxes in eclipse computed from Kappa fit at all potentials. Integral extends from 1 eV to 
100 keV. 


Sept 01; Eclipse; Kappa fits 
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Figure 122. Ion fluxes in eclipse computed from Kappa fit at all potentials. Integral extends from 1 eV to 100 
keV. 
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Sept 01; Eclipse; Kappa & Kappa fits 



Fluxes (1/(cm 2 s sr)) 


Figure 123. Net fluxes in eclipse computed from Kappa fits at all potentials. Integral extends from 1 eV to 100 
keV. 


Sept 01; Eclipse; Kappa & Kappa fits 
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Figure 124. Measured potential as a function of the potential computed from the minimum in the net flux, 
where the fluxes are computed from Kappa fits to the electron and ion fluxes. 
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Sept 01 ; Eclipse; Max & Single Max fits 
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Figure 125. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Maxwellian fit to the electron flux and a Single Maxwellian fit to the ion flux. 


Sept 01 ; Eclipse; Max & Double Max fits 
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Figure 126. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Maxwellian fit to the electron flux and a Double Maxwellian fit to the ion flux. 
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Sept 01 ; Eclipse; Max & Kappa fits 
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Figure 127. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Maxwellian fit to the electron flux and a Kappa fit to the ion flux. 


Sept 01 ; Eclipse; Single Max & Max fits 
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Figure 128. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Single Maxwellian fit to the electron flux and a Maxwellian fit to the ion flux. 
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Sept 01 ; Eclipse; Single Max & Double Max fits 
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Figure 129. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Single Maxwellian fit to the electron flux and a Double Maxwellian fit to the ion flux. 


Sept 01 ; Eclipse; Single Max & Kapp fits 
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Figure 130. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Single Maxwellian fit to the electron flux and a Kappa fit to the ion flux. 
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Sept 01 ; Eclipse; Double Max & Max fits 
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Figure 131. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Double Maxwellian fit to the electron flux and a Maxwellian fit to the ion flux. 


Sept 01 ; Eclipse; Double Max & Single Max fits 
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Figure 132. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Double Maxwellian fit to the electron flux and a Single Maxwellian fit to the ion flux. 
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Sept 01; Eclipse; Double Max & Kappa fits 
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Figure 133. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Double Maxwellian fit to the electron flux and a Kappa fit to the ion flux. 


Sept 01; Eclipse; Kappa & Max fits 
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Figure 134. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Kappa fit to the electron flux and a Maxwellian fit to the ion flux. 
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Sept 01; Eclipse; Kappa & Single Max fits 



Sept 01 ; Eclipse; Kappa & Single IV 


548.2 


> 

ro 

c 

0) 

o 

Q. 


141.3 


37.6 


3.59 



0 25 50 

Number no solution 


Figure 135. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Kappa fit to the electron flux and a Single Maxwellian fit to the ion flux. 


Sept 01 ; Eclipse; Kappa & Double Max fits 
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Figure 136. Measured potential computed from the minimum in the net flux, where the fluxes are computed 
from a Kappa fit to the electron flux and a Double Maxwellian fit to the ion flux. 
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Sept 01; Eclipse; Kappa & Max fits 
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Figure 137. Net fluxes in eclipse at all potentials with electron fluxes computed from Kappa fit and ion fluxes 
computed from Kappa fit. Integral extends from 1 eV to 100 keV. 
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Table 2. Quality of potential predictions made from various fits to measured incident fluxes. 


Electron fit 

Ion fit 

Number 
with no 
match 

Average 

error 

Number Number 

within factor within factor 
of 1.5 of 3 

Maxwellian 

Maxwellian 

338 

0.37 

384 

599 

Maxwellian 

Single Maxwellian 

182 

0.897 

288 

604 

Maxwellian 

Double Maxwellian 

307 

0.479 

322 

527 

Maxwellian 

Kappa 

263 

0.478 

187 

344 

Single Maxwellian 

Maxwellian 

194 

0.669 

564 

692 

Single Maxwellian 

Single Maxwellian 

152 

1.466 

449 

684 

Single Maxwellian 

Double Maxwellian 

243 

0.777 

349 

605 

Single Maxwellian 

Kappa 

304 

0.1417 

169 

343 

Double Maxwellian 

Maxwellian 

239 

0.642 

484 

654 

Double Maxwellian 

Single Maxwellian 

207 

1.22 

346 

70 

Double Maxwellian 

Double Maxwellian 

257 

0.707 

300 

597 

Double Maxwellian 

Kappa 

251 

0.1535 

182 

375 

Kappa 

Maxwellian 

102 

0.362 

634 

774 

Kappa 

Single Maxwellian 

69 

0.388 

455 

822 

Kappa 

Double Maxwellian 

134 

0.36 

459 

712 

Kappa 

Kappa 

165 

0.112 

230 

457 
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Table 3. Comparison of the accuracy of potential predictions from the fits when the flux integrals 
was done using only those energy bins used to create the flux and when the flux integrals 
extended from 1 eV to 100 keV. 


Electron fit 

Ion fit 

Number with 

Bins summed no match 

Average 

error 

Number 
within factor 
of 1.5 

Maxwellian 

Maxwellian 

1 to 100 keV 

338 

0.37 

384 

Maxwellian 

Maxwellian 

Same as data 

342 

0.84 

336 

Single Maxwellian 

Single Maxwellian 

1 to 100 keV 

152 

1.466 

449 

Single Maxwellian 

Single Maxwellian 

Same as data 

152 

1.71 

393 

Double Maxwellian 

Double Maxwellian 

1 to 100 keV 

257 

0.707 

300 

Double Maxwellian 

Double Maxwellian 

Same as data 

133 

0.778 

494 

Kappa 

Maxwellian 

1 to 100 keV 

102 

0.362 

634 

Kappa 

Maxwellian 

Same as data 

106 

0.469 

665 

Kappa 

Kappa 

1 to 100 keV 

165 

0.112 

230 

Kappa 

Kappa 

Same as data 

70 

0.177 

298 
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